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SECTION  1.  INTRODUCTION 


The  Open  Burn/Open  Detonation  Dispersion  Model  (OBODM)  is  intended  for 
use  in  evaluating  the  potential  air  quality  impacts  of  the  open-air  burning  or 
detonation  of  obsolete  munitions,  solid  rocket  propellants,  and  manufacturing 
wastes  at  Department  of  Defense  and  Department  of  Energy  installations.  OBODM 
predicts  the  downwind  transport  and  dispersion  of  pollutants  using  cloud  rise  and 
dispersion  model  algorithms  taken  from  existing  dispersion  models.  OBODM  is 
strictly  a  transport  and  dispersion  model;  it  does  not  contain  a  source  preprocessor 
to  predict  the  total  quantities  of  pollutants  released  by  an  open  burn  or  detonation. 
Rather,  it  is  designed  to  use  either  empirical  emissions  factors  such  as  those 
derived  from  measurements  in  the  Dugway  Proving  Ground  Bang  Box™  or  the 
emissions  predicted  by  a  products  of  combustion  model  (for  example,  Baroody, 
1994).  (OBODM  contains  a  data  base  of  empirical  BangBox  emissions  factors.) 
OBODM  calculates  peak  concentration,  time-mean  concentration,  and  dosage  (time- 
integrated  concentration)  for  quasi-continuous  (open  burn)  and  instantaneous  (open 
detonation)  releases.  It  can  also  consider  the  effects  on  concentration  and  dosage 
of  the  gravitational  settling  and  deposition  of  particulates  such  as  soil  particles 
entrained  into  the  cloud  from  an  open  detonation. 

OBODM  can  consider  two  types  of  quasi-continuous  or  instantaneous 
sources:  volume  and  line.  (The  model  automatically  represents  a  line  source  by  a 
series  of  volume  sources.)  The  vast  majority  of  open  burn/open  detonation 
(OB/OD)  releases  can  be  modeled  using  volume  sources.  Line  sources  are  available 
to  model  relatively  unusual  OB/OD  source  configurations  or  other  types  of  releases 
such  as  dust  released  from  the  detonations  used  in  surface  mining  operations. 

Volume  I  of  this  two-volume  OBODM  User's  Guide  provides  user's  instruc¬ 
tions,  including  installation  procedures,  a  tutorial  on  model  use,  examples  of  all 
menus,  and  example  problems.  This  volume  provides  a  detailed  technical 
description  of  OBODM,  including  the  mathematical  models  and  the  procedures  used 
to  parameterize  meteorological  conditions  in  the  surface  mixing  layer. 


1 


Blank  Page 


SECTION  2.  MATHEMATICAL  MODELS 


2.1  PLUME  AND  CLOUD  RISE 

2.1.1  Plume  Rise  (Open  Burn) 

OBODM  uses  the  Briggs  (1971)  plume  rise  equations  to  predict  the  buoyant 
rise  of  the  plume  from  an  open  burn.  Following  the  derivation  given  by  Dumbauld 
et  al.  (1973),  these  equations  have  been  modified  to  account  for  the  horizontal 
dimensions  of  a  very  large  burn  (diameter  greater  than  about  5  m).  With  a  stable 
thermal  stratification,  the  buoyant  rise  Ahcs(x)  at  downwind  distance  x  from  a  burn 
of  greater  than  15  s  duration  is  given  by 


(  r-  P 

3  F 


Ahc,(x)  = 


uy;S 


[1  -  cos  (S'12  x/u)]  + 


f  > 

3' 

1/3 

(■  \ 

_ 

rjL\-x<mS-m 

\yj 

vyj 

6F 
uy2  S 


“1 1/3 


\YCJ 


;  x  >  7vuS 


-1/2 


(2-1) 


gHeR 

ft  P p  T  a 


(2-2) 


g_  50 
Ta  dz 


(2-3) 


where 

u  =  mean  wind  speed  (m  s'1)  at  release  height  h  or  2  m,  whichever  is 
greater 

yc  =  continuous  source  entrainment  coefficient  (default  value  is  0.6) 
rR  =  radius  of  the  burning  material  (m) 
g  =  acceleration  due  to  gravity  (9.8  m  s'2) 

Hc  =  material  effective  heat  content  (cal  g'1) 
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R  =  material  burn  rate  (g  s'1) 

cp  =  specific  heat  of  air  at  constant  pressure  (0.24  cal  g'1  K !) 

Pa  =  air  density  (g  m‘3) 

Ta  =  air  temperature  (K) 

^  =  vertical  potential  temperature  gradient  in  layer  through  which  cloud 
rises  (K  m'1) 


With  an  unstable  or  neutral  thermal  stratification  (S6/3Z  ^  0),  buoyant  rise  at 
downwind  distance  x  from  a  burn  of  greater  than  15s  duration  is  given  by 
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Inspection  of  Equations  (2-1)  through  (2-3)  shows  that  the  predicted  stable  rise 
approaches  infinity  in  the  limit  as  the  vertical  potential  temperature  gradient  S0/3z 
approaches  zero.  Thus,  it  is  possible  for  the  plume  rise  predicted  by  Equation  (2-1) 
to  exceed  that  predicted  for  adiabatic  or  unstable  conditions  by  Equation  (2-4). 
Consequently,  OBODM  uses  the  smaller  of  the  plume  rises  predicted  by  Equations 
(2-1)  and  (2-4)  under  stable  conditions. 

2.1.2  Cloud  Rise  (Open  Detonation) 

OBODM  calculates  the  buoyant  rise  of  the  cloud  from  a  detonation  or  burn  of 
duration  less  than  1  5  s  using  equations  from  the  Rocket  Exhaust  Effluent  Diffusion 
Model  (REEDM)  (Bjorklund  et  al.,  1982).  These  equations  were  derived  by 
Dumbauld  et  al.  (1973)  following  the  same  reasoning  as  used  by  Briggs  (1971)  in 
the  derivation  of  his  rise  equations  for  a  continuous  source.  The  buoyant  cloud  rise 
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Ah,s(x)  at  downwind  distance  x  with  a  stable  stratification  is  given  by 
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where  Wd  is  the  quantjtV  of  material  detonated  in  grams,  Yi  is  the  instantaneous 
source  entrainment  coefficient,  and  the  other  terms  are  as  defined  above.  Based  on 
the  recommendations  of  Bjorklund  et  al.  (1982),  the  default  value  for  Yi  is  0.64. 
Although  Dumbauld  et  al.  (1973)  derived  an  equation  analogous  to  Equation  (2-4) 
for  an  instantaneous  source,  Bjorklund  et  al.  (1982)  use  the  stable  rise  equation  of 
Equation  (2-7)  to  approximate  instantaneous  source  cloud  rise  with  adiabatic  or 
unstable  conditions  by  setting  the  potential  temperature  gradient  equal  to  3.344  x 
10~4  K  m'1.  This  potential  temperature  gradient  was  empirically  determined  using 
observations  of  the  buoyant  rise  of  the  exhaust  clouds  from  large  solid  rocket 
motors  such  as  those  used  by  the  Space  Shuttle  and  Titan  34D.  Thus,  it  is  likely 
that  this  empirical  value  implicitly  includes  the  buoyancy  flux  Fi  of  these  large 
sources. 

As  an  alternative  to  the  use  of  Equation  (2-7)  when  there  is  an  adiabatic  or 
unstable  lapse  rate,  OBODM  allows  the  user  to  select  the  adiabatic  cloud  rise 
equation  derived  by  Dumbauld  et  al.  (1973).  If  the  adiabatic  rise  equation  is  used, 
the  downwind  distance  to  final  cloud  rise  is  defined  in  a  manner  similar  to  the 
distance  to  final  plume  rise  for  a  continuous  source  [see  Equations  (2-5)  and  (2-6)]. 
The  resulting  equations  are 
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Note  that  the  xmax  given  by  Equation  (2-1 1)  is  analogous  to  the  3.5  x*  given  by 
Equations  (2-5)  and  (2-6).  The  derivation  of  Equation  (2-1 1)  is  given  in  Appendix 
A. 


It  should  be  noted  that  the  cloud  from  an  open  burn  does  not  always  behave 
as  a  quasi-continuous  source.  Specifically,  if  the  burn  time  is  very  short,  the 
resulting  cloud  rises  and  disperses  in  the  same  manner  as  an  instantaneous  source. 
The  burn  time  tB  is  given  by 
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where  Wb  is  the  quantity  burned  in  grams.  If  tB  is  less  than  or  equal  to  15  s, 
OBODM  treats  the  source  as  an  instantaneous  release  with  Wb  substituted  for  Wd 
in  Equation  (2-8). 


2.1.3  Inversion  Penetration 


OBODM  accounts  for  the  fact  that  a  plume  or  cloud  may  partially  or 
completely  penetrate  the  elevated  layer  above  the  top  of  the  surface  mixing  layer 
using  a  technique  similar  to  that  used  by  the  PPSP  model  (Weil  and  Brower,  1  984). 
Although  this  technique  was  specifically  developed  to  address  plume  penetration 
into  the  elevated  stable  layer  above  a  convective  or  adiabatic  surface  mixing  layer, 
it  is  also  applied  by  OBODM  to  the  case  of  the  layer  capping  a  stable  surface  mixing 
layer. 

Plume  or  cloud  height  H(x)  at  downwind  distance  x  is  given  by 

H(x)  =  h  +  A  h(x)  (2-13) 

where  h  is  the  release  height  above  ground  level  (normally  zero)  and  Ah(x)  is  the 
plume  rise  obtained  from  Equation  (2-1),  (2-4),  (2-7),  or  (2-9),  depending  on  the 
source  type  and  stability.  If  Ah  exceeds  0.67(Hm-h),  where  Hm  is  the  mixing  depth, 
a  second  plume  rise  Ah'(x)  is  calculated  using  Equation  (2-1)  or  (2-7),  depending  on 
source  type,  and  a  vertical  potential  temperature  gradient  dQld z  of  0.01  K  m1.  (If 
the  potential  temperature  gradient  used  in  the  original  rise  calculation  exceeds  0.01 
K  rrf1,  Ah'(x)  is  set  equal  to  the  original  Ah(x).)  The  fraction  P  of  the  plume  or  cloud 
in  the  elevated  layer  is  then  given  by 


6 


1 


p  = 


(Hm-h 

Ah'(x) 


<  0.5 


1.5 

0 


„  h  (Hm-h 

m  ;  0.5  <  ■  m.  <  15 


A  h\x) 


;  1.5  < 


Ah'  (x) 
A  h'{x) 


(2-14) 


OBODM  assumes  that  the  portion  of  a  cloud  or  plume  that  stabilizes  above 
the  top  of  the  surface  mixing  layer  does  not  contribute  to  concentrations  within  the 
mixing  layer.  The  effective  source  strength  Qe(x)  for  the  plume  or  cloud  material 
within  the  mixing  layer  is  defined  by 

Q,(x)  =  Q0(i-P)  (2'15) 

ii 

where  Q0  is  the  original  source  strength  (quantity  of  pollutant  released  by  the  burn 
or  detonation).  The  effective  plume  or  cloud  rise  for  the  material  within  the  mixing 
layer,  which  is  substituted  for  Ah(x)  in  Equation  (2-13),  is  defined  as 

Ahe(x)  =  (0.62  +  0.38 P)  (H,„-h)  (2-16) 


2.1 .4  Plume/Cloud  Rise  for  Line  Sources 


OBODM  is  not  specifically  designed  to  calculate  buoyant  plume  or  cloud  rise 
for  a  line  source.  The  model  can  be  directed  to  use  the  volume  source  rise 
equations  described  above  to  estimate  buoyant  rise  for  a  line  source  by  entering  a 
nonzero  heat  content  Hc,  but  the  user  should  recognize  that  the  resulting  rise 
predictions  may  not  be  realistic. 

2.2  PARTICLE  DISTRIBUTION,  SETTLING  VELOCITY,  AND  REFLECTION 
COEFFICIENT 

OBODM  requires  the  particle  size  distribution  in  order  to  consider  the  effects 
of  gravitational  settling  on  concentration,  dosage,  and  dry  deposition.  The  size 
distribution  is  divided  into  N  size  categories.  For  each  category,  the  model  requires 
the  diameter,  settling  velocity,  fraction  of  total  mass  contained  in  the  category,  and 
surface  reflection  coefficient.  The  user  may  enter  all  of  this  information  or  direct 
OBODM  to  compute  it  from  the  mass  median  particle  diameter,  geometric  standard 
deviation  of  the  particle  size  distribution,  and  density  of  the  particulates. 

2.2.1  Particle  Size  Distribution 


The  particle  size  distribution  automatically  generated  by  OBODM  is  a  log¬ 
normal  distribution  that  extends  2.3  standard  deviations  to  each  side  of  the  median. 
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The  particle  diameter  at  the  lower  boundary  of  the  nth  size  category  is  given  by 

In  ad 


dm  exp 


dn 


4. 6(n  - 1) 


N 


<?d3 


(2-17) 


where 

dm  =  mass-median  diameter  (m) 

ad  =  geometric  standard  deviation  (od  >  1) 

N  =  number  of  particle  size  categories  (maximum  of  20) 

t 

if 

If  Od  is  entered  as  less  than  or  equal  to  unity,  OBODM  assumes  monodisperse 
particles  of  diameter  dm. 

The  mass  fraction  of  particles  in  the  nth  size  category  is  given  by 


h- 


0.01  +-A„ 


1 


n  =  1  or  N 


> 


1  <  n  <  N 


(2-18) 


where 


erf 


In  d„+i  -  In  d„ 


42  a. 


■erf 


^ln  d„  -  In  d  ^ 


42  G, 


(2-19) 


Whether  the  particle  size  distribution  is  automatically  generated  from  a  log¬ 
normal  distribution  or  entered  by  the  user,  OBODM  represents  each  particle  size 
category  by  its  mass-mean  diameter,  which  is  given  by 


dn 


dn+i  "l-  dndn+i  dndn+i  d a 

4 


-y/3 


(2-20) 
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2.2.2  Gravitational  Settling  Velocity 


OBODM  uses  the  method  of  McDonald  (1960)  to  compute  the  terminal 
settling  velocity  for  the  mass-mean  diameter  of  each  particle  size  category.  The 
terminal  fall  velocity  of  the  nth  size  category  is  given  by 


v  = 

r  sn 


pg(dny 


18  Pa 


;  Q"  <24  ( Stokes  law) 


Pa 


Pa 


Ri 


91 

Q'i 


;  24  <  Q"  ;  Q't  <  Q"  <  Q'M 


(2-21) 


Q"  = 


400gppa  ( dj 3 

3K 


(2-22) 


In 


<1  = 


Rj+i 
Ri ) 


In 


v  Qi  > 


(2-23) 


where  the  mass-mean  particle  dn  is  in  centimeters,  pa'  is  the  air  density  in  grams  per 
cubic  centimeter  (pa  is  used  in  Section  2.1  to  indicate  air  density  in  grams  per  cubic 
meter),  and 

p  =  particle  density  (g  cm'3) 
g  =  acceleration  due  to  gravity  (m  s'2) 

R|,  Qj'  =  coefficients  listed  in  Table  1 

pa  =  absolute  viscosity  of  air  (g  cm'1  s'1) 

=  7-6342 Xj^  f_T^y/2  24  (2-24) 

(Ta  +  120)  \  296.16) 

Ta  =  air  temperature  (K) 

The  logarithmic  interpolation  shown  in  Equation  (2-21)  for  Q"  greater  than  24 
depends  on  an  empirical  relationship  between  Q"  and  the  Reynolds  number  Re. 

After  Q"  has  been  calculated,  the  values  of  Q'i  and  Q'i+i  bounding  Q"  and  the 
corresponding  values  of  Ri  and  Ri+i  are  selected  from  Table  1. 
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Table  1.  Values  of  Ri  and  Q'i  for  Settling  Velocity  Computations. 

i  Ri  Q'i 


1 

1 

24 

2 

2 

64 

3 

4 

140 

4 

6 

230 

5 

10 

420 

6 

20 

1 .1 6x1 03 

7 

r  1 

40 

2.88x1 03 

8  V 

60 

5.40x1 03 

9 

100 

1 .20x1 04 

10 

200 

3.20x1 04 

1 1 

400 

9.76x1 04 

12 

600 

1 .94x1 05 

13 

IxlO3 

4.60x1 05 

14 

2x1 03 

1 .64x1 06 

15 

4x1 03 

6.40x1 06 

16 

6x1 03 

1 .44x1 07 

17 

IxlO4 

4. 1 0x1 07 

18 

2x1 04 

1 .76x1 08 

19 

4x1 04 

7.33x1 08 

20 

6x1 04 

1 .65x1 09 

21 

IxlO5 

4.40x1 09 

22 

2x1 05 

1 .64x1 010 
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2.2.3  Reflection  Coefficient 


The  reflection  coefficient  is  the  fraction  of  material  reaching  the  surface 
through  the  combined  processes  of  gravitational  settling  and  atmospheric 
turbulence  that  is  assumed  to  be  reflected  from  the  surface.  It  ranges  from  unity 
for  complete  reflection  to  zero  for  complete  deposition.  Thus,  an  upper  bound  on 
dosage  or  concentration  can  be  obtained  by  setting  the  reflection  coefficient  equal 
to  unity,  while  the  effects  of  depletion  by  dry  deposition  can  be  maximized  by 
setting  the  reflection  coefficient  to  zero.  If  the  user  does  not  enter  reflection 
coefficients,  OBODM  automatically  estimates  them  using  a  semi-empirical 
relationship  (Dumbauld  et  al.,  1976)  between  the  reflection  coefficient  and  settling 
velocity.  Rather  than  require  the  user  to  estimate  the  reflection  coefficient  from  the 
Dumbauld  et  al.  (1976)  nomogram,  OBODM  uses  the  Bjorklund  et  al.  (1982)  fit  to 
the  nomogram 


i/ 


O  ;  V„  >03 ms~' 

2.5(03  -  VJ  ;  0.3  <VW  <0.05  ms~'  ■ 
1.0-1.0179  V'J;  ;  K„,  <0.05 


(2-25) 


Note  that  OBODM  requires  a  surface  reflection  coefficient  that  does  not  vary  with 
downwind  distance.  The  Overcamp  (1976)  reflection  coefficient,  which  is  a 
function  of  deposition  velocity  and  downwind  distance,  is  not  suitable  for  use  in 
OBODM. 


2.3  CONCENTRATION 


OBODM  uses  the  concentration  algorithms  from  Dugway  Proving  Ground's 
RTVSM  model  (Bjorklund,  1990).  These  algorithms  are  based  on  the  same  steady- 
state  Gaussian  dispersion  model  concepts  used  by  many  other  models.  OBODM 
has  two  basic  source  types,  the  Gaussian  puff  and  a  square-wave  quasi-continuous 
source.  The  square-wave  quasi-continuous  source  model  can  be  derived  by  inte¬ 
grating  the  Gaussian  puff  equation  with  respect  to  time  over  the  limits  0  to  x, 
where  T  is  the  release  duration.  Thus,  the  quasi-continuous  source  model  is  the 
analytic  equivalent  of  representing  a  quasi-continuous  release  by  an  overlapping 
series  of  Gaussian  puffs. 


The  concentration  produced  by  an  instantaneous  or  quasi-continuous  volume 
source  at  downwind  distance  x,  crosswind  distance  y,  height  z  above  the  surface, 
and  time  t  after  source  function  is  given  by 

X(x,y,z,t)  =  Xp(x)  •  VT(x,z)  •  LT(x,y)  •  AT(x,t)  •  DC(x)  (2-26) 


li 


where  Xp,  VT,  LT,  AT,  and  DC  are  the  peak  concentration,  vertical,  lateral, 
alongwind,  and  decay  terms,  respectively. 


The  peak  concentration  term  is  defined  as 

K  Qe  (X) 


ZP  (X)  = 


(2x) 


3/2 


0-xCTv<7: 


K  Qe  {X) 


(2n)  cr  cr.  u  r 


;  instantaneous  source 


;  quasi  -  continuous  source 


(2-27) 


where 

K  =  units  conversion  factor 

r 

::j 

Qe(x)  =  effective  source  strength  [see  Equation  (2-15)] 

ax  =  alongwind  (longitudinal)  dispersion  coefficient  (standard 
deviation  of  alongwind  concentration  distribution)  (m) 

av  =  lateral  dispersion  coefficient  (standard  deviation  of  crosswind 
concentration  distribution)  (m) 

a2  =  vertical  dispersion  coefficient  (standard  deviation  of  vertical 
concentration  distribution)  (m) 

u  =  mean  transport  wind  speed  (m  s'1) 

T  =  dissemination  duration  (burn  time  tB)  (s) 


The  transport  wind  speed  and  dispersion  coefficients  are  discussed  in  Section  2.6. 

The  vertical  term  for  gases  and  aerosols  with  negligible  settling  velocities  is 
given  by 


VT  (x,z)  =  £ 


i=0 


exp 


2 


f2iH,„  -  H+zV 


+ 


exp 


r2iH„+  H  +  zV 
V  cr. 


(2-28) 


+  I 


i=l 


exp 


f2  iHm  +H-zV' 


cr.  J 


+  exp 


( 2iHm-  H-z}1 
\  cr. 
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where  Hm  is  the  mixing  depth  and  H  is  the  effective  release  height  [see  Equations 
(2-1 3)  and  (2-1 6)].  At  the  downwind  distance  where  the  exponential  terms  for  i 
equals  3  exceed  exp(-IO),  the  cloud  or  plume  is  approximately  uniformly  mixed 
within  the  surface  mixing  layer.  Equation  (2-28)  is  then  replaced  by 


VT(x,z)  = 


(2-29) 


which  results  in  a  box  model  with  uniform  vertical  mixing. 

OBODM  can  also  consider  the  effects  on  concentration  and  dosage  of  the 
gravitational  settling  and  dry  deposition  of  particulates  with  appreciable 
gravitational  settling  velocities  by  using  a  different  and  more  generalized  vertical 
term.  The  derivations  of  this  algorithm  and  the  associated  surface  deposition 
algorithms  are  discussed  by  Cramer  and  Bowers  (1980).  Similar  algorithms  are  also 
used  in  other  models  such  as  the  original  ISC  (Bowers  et  al.,  1979),  a  regulatory 
model  designed  for  application  to  industrial  particulate  emissions,  and  FSCBG 
(Teske  et  al.,  1993),  an  aerial  spray  dispersion  and  deposition  model. 


As  discussed  in  Section  2.2,  OBODM  allows  the  user  to  divide  the  source 
strength  for  particulates  with  appreciable  gravitational  settling  velocities  into  as 
many  as  20  settling-velocity  categories,  or  the  model  can  be  directed  to  generate  a 
log-normal  distribution  from  the  mass-median  diameter  and  geometric  standard 
deviation.  For  particulates  in  the  nth  settling-velocity  category,  the  vertical  term  is 


r. 

r 

VT  (x,z)=  ^  • 

£ 

4  exp 

/  =  0 

2  V 

2/  Hm  -  H  +  z  +  Vw  x  /  u 
o'.  , 


+  4+l/  exP 


1  f  2i  Hm  +  H  +  z  -  V  x/u  V 


21 


cr 


+  1 


4  exp 


]_ 

2 


2  i  Hm  +  H  -  z  —  Vm  x  /  u 
<  cr;  , 


(2-30) 


+  4'1  exp 


1  f  2  i  Hm  -  H  -  z  +  Vm  x/u'2 


V 


cr. 


J 


13 


where 

^  =  mass  fraction  of  source  strength  in  nth  settling-velocity  category 

Vsn  =  settling  velocity  of  particulates  in  the  n,h  settling-velocity 
category  (m  s'1)  (see  Section  2.2.2) 

Sn  =  surface  reflection  coefficient  for  particulates  in  the  nth 
settling-velocity  category  (see  Section  2.2.3) 

For  convenience,  0°  is  defined  to  be  unity  in  Equation  (2-30).  The  total 
concentration  is  obtained  by  summing  over  the  N  settling-velocity  categories. 

OBODM  is  a  surface  depletion  deposition  model  that  selectively  depletes 
material  at  the  ground  rather  than  through  the  entire  vertical  extent  of  the  cloud. 
Because  of  the  selective*  depletion  at  the  surface,  uniform  vertical  mixing  is  never 
achieved  and  Equation  (2-30)  cannot  be  replaced  with  the  box  model  of  Equation 
(2-29). 


The  OBODM  lateral  term  LT  and  decay  term  DC  are  respectively  given  by 


LT(x,y)  =  exp 


— 

f  \ 

2 

1 

y 

'  2 

\  (Xy/ 

(2-31) 


and 


DC(x)  =  exp  [-kx/  u] 


(2-32) 


where  k  is  the  decay  coefficient  (fraction  of  the  material  lost  per  unit  time).  For 
example,  if  T 1/2  is  the  pollutant  half  life  in  seconds,  k  can  be  obtained  from 


0.693 


Tin 


(2-33) 


The  alongwind  term  is  given  by 


AT(x,t)  = 


exp 


fx-ut^' 


vo-,; 


;  instantaneous  source 


(2-34) 


erf 


{  x  -  u(t- 


O',  J 


-erf 


^  x-ut^ 


\42<jxJ 


;  quasi  -  continuous  source 
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where  erf  is  the  Gaussian  error  function.  If  OBODM  is  used  to  calculate 
concentration  time  histories  at  selected  (discrete)  receptors,  concentrations  are 
calculated  at  each  receptor  at  equal  time  steps  between  times  ti  and  t2,  which  are 
defined  as 


x  2.45  <7X 
u 


(2-35) 


^2 


*  +  2.45  £7 

- = - -  ;  instantaneous  source 

u 

x  +  2.45  a 

x  + - _ - -  ;  quasi  -  continuous  source 

u 


(2-36) 


In  addition  to  calculating  the  concentration  time  histories  at  discrete  receptors, 
OBODM  can  be  directed  to  calculate  peak  or  maximum  time-mean  concentrations  at 
all  receptors.  If  peak  concentrations  are  selected,  the  alongwind  term  of  Equation 
(2-34)  is  evaluated  at  time  tp,  defined  as 


t 


p 


x 

—  ;  instantaneous  source 

u 


X  X 

—  +  —  ;  quasi  -  continuous  source 
.  u  2 


(2-37) 


If  time-mean  concentration  is  specified,  the  alongwind  term  for  an  instanta¬ 
neous  source  is  given  by 


AT(x)  = 


f 

V 


u  xA 


\  ( 


rf{ 


u  X  A 

2^2  ax  a 


(2-38) 


where  ta  is  the  concentration  averaging  time.  The  equation  corresponding  to 
Equation  (2-38)  for  a  quasi-continuous  source  is 
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AT(x)  = 


4i  u 


ta  J 


u(r-fA)  /  u  (j-T/S 
4l  crx  V  2V2  <jx  j 


u(r  +  rA)  erf  f  uQ  +  Ta)^ 
V2  crx  I  2\/2  crx  > 


2 


f 

(  u(t-ta) 

A 

exp 

\ 

V  2V2  (7X/ 

y 

(2-39) 


f 


-  exp 


u(^  +  ta) 
2V2  crxy 


x/  J 


Note  that  the  use  of  Equation  (2-38)  or  (2-39)  results  in  a  concentration  at  each 
receptor  that  is  the  mean  concentration  over  an  interval  of  length  xA  centered  on  the 
time  of  the  peak  concentration  at  that  receptor.  These  equations  can  be  derived  by 
integrating  the  generalized  concentration  equation  with  respect  to  time  over  the 
limits  of  tp-xA/2  to  tp  +  xA/2  and  dividing  by  xA. 

2.4  DOSAGE 

2.4.1  Instantaneous  Source  (Open  Detonation) 


The  OBODM  dosage  algorithms  are  also  based  on  conventional  steady-state 
Gaussian  dispersion  model  concepts.  The  dosage  (time-integrated  concentration) 
produced  by  an  instantaneous  volume  source  at  downwind  distance  x,  crosswind 
distance  y,  and  height  z  above  the  surface  is  given  by 

D(x,y,z)  =  Dp(x)  •  VT(x,z)  •  LT(x,z)  •  DC(x)  (2-40) 


where  Dp,  VT,  LT,  and  DC  are  the  peak  dosage,  vertical,  lateral,  and  decay  terms, 
respectively.  The  equations  for  the  vertical,  lateral,  and  decay  terms  are  identical 
to  the  equations  previously  given  in  Section  2.3.  The  peak  dosage  term  is 


DP(x)  = 


K  Q,(x)  _ 
2n  <jy  <7Z  u 


(2-41) 


where  the  terms  are  also  as  defined  in  Section  2.3. 
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2.4.2  Quasi-Continuous  Source  (Open  Burn) 


OBODM  computes  the  dosage  produced  by  a  quasi-continuous  source  by 
numerical  integration  of 


‘2  ..  '  ‘  ‘ 

D(x,y,z)  =  |  x(x,y,z,t)dt'  (2-42) 

‘1 


where  the  concentration  %  is  given  by  Equation  (2-26)  and  the  limits  of  integration 
ti  and  t2  are  given  by  Equations  (2-35)  and  (2-36),  respectively. 

2.5  DEPOSITION 

2.5.1  Gravitational  Deposition 


Particulates  in  OB/OD  emissions  can  consist  of  those  contained  in  the 
products  of  combustion  (for  example,  aluminum  oxide  particulates  for  some  types 
of  solid  rocket  propellants),  soil  dust  entrained  by  detonations,  and  fragments  of 
munition  casings.  OBODM  predicts  the  gravitational  settling  and  deposition  (i.e., 
fallout)  of  particulates  with  appreciable  settling  velocities  using  an  algorithm  of  the 
same  type  as  used  by  several  models  including  REEDM,  RTVSM,  and  the  original 
ISC.  The  vertical  term  of  Equation  (2-30)  accounts  for  the  effects  on  ambient 
particulate  concentrations  of  gravitational  settling  and  deposition  by:  (1)  depleting 
the  airborne  mass  to  account  for  the  losses  at  the  surface,  and  (2)  changing  the 
shape  of  the  vertical  concentration  distribution.  If  the  settling  velocities  of  the 
particulates  are  sufficiently  small  that  a  fraction  of  the  total  mass  is  initially  mixed 
upward  to  the  top  of  the  surface  mixing  layer  where  it  is  reflected  downward, 
OBODM  predicts  a  secondary  maximum  in  the  concentration  profile  as  this  fraction 
reaches  the  surface.  This  phenomena  can  be  observed  in  field  experiments  with 
sufficiently  small  particulates  if  the  measurements  extend  far  enough  downwind 
(see  Dumbauld  et  al.,  1976). 

The  OBODM  gravitational  deposition  algorithm  can  be  derived  from  Equation 
(2-26)  with  the  vertical  term  of  Equation  (2-30)  (see  Cramer  and  Bowers,  1980). 
For  both  instantaneous  and  quasi-continuous  sources,  deposition  at  downwind 
distance  x  and  crosswind  distance  y  for  particulates  in  the  nth  settling  velocity 
category  is  given  by 


DEPn (x, y)  =  DEP„p(x)  •  VTnD  •  LT(x)  •  DC(x)  (2-43) 


where  DEPnp,  VTnD,  LT,  and  DC  are  the  peak  deposition,  vertical  (deposition),  lateral, 
and  decay  terms,  respectively.  The  lateral  and  decay  terms  are  as  defined  in 
Section  2.3.  The  peak  deposition  term  for  the  nth  settling-velocity  category  is  given 
by 
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(2-44) 


DEP„p(x) 


K  Qe(x)(l-8„) 
2tt  (Jy  <7;  X" 


x"  =  X  +  Xz  -  Xrz(l-P) 


(2-45) 


where  xz  is  the  vertical  virtual  distance,  xrz  is  the  rectilinear  vertical  expansion 
distance,  and  P  is  the  vertical  expansion  coefficient.  These  three  parameters, 
which  are  used  to  calculate  the  vertical  dispersion  coefficient  Gz,  are  discussed  in 
Section  2.6.  Note  that,  except  for  a  source  with  a  very  large  initial  vertical 
dimension,  the  x"  given  by  Equation  (2-45)  is  approximately  equal  to  the  downwind 
distance  x. 

t  ' 

The  deposition  vertical  term  for  particulates  in  the  n,h  settling-  velocity 
category  is  given  by 


VTllD(x)  =  </>„  { M  +  N} 


(2-46) 


where 


M  = 


exp 


H  -  V  m  x  /  u 
V  cr_-  J 


V 


(2-47) 
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1_ - 
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V  2  v 


j3(2iHm  +  H)  + 
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v  .  x" 


J 


) 


exp 
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(2-48) 
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As  in  the  case  of  Equation  (2-30),  0°  is  defined  as  unity  in  Equation  (2-48).  OBODM 
obtains  the  total  gravitational  deposition  by  summing  over  the  N  settling-velocity 
categories. 

2.5.2  Dry  Deposition 

Gravitational  deposition  is  not  the  only  mechanism  by  which  material  from  an 
OB/OD  cloud  or  plume  can  be  transferred  to  the  surface.  Small  quantities  of  gases 
and  aerosols  that  come  in  contact  with  the  surface  can  be  retained  as  a  result  of 
processes  such  as  absorption,  vegetative  uptake,  and  impaction  on  surface 
roughness  elements.  This  overall  process  is  usually  described  as  dry  deposition. 

For  a  short-term  release,  the  total  dry  deposition  at  downwind  distance  x  and 
crosswind  distance  y  is  typically  estimated  from 

DRY(x,y)  =  vd  •  D(x,y,0)  (2‘49) 

where  D(x,y,0)  is  the  total  ground-level  dosage  at  that  point  and  vd  is  an  empirical 
or  theoretical  dry  deposition  velocity.  A  summary  of  empirical  deposition  velocities 
for  gases  and  small  particles  is  given  by  McMahon  and  Denison  (1982). 

Gaussian  dispersion  models  that  use  the  deposition  velocity  parameterization 
normally  use  an  effective  source  strength  that  decreases  with  downwind  distance 
to  account  for  the  losses  due  to  dry  deposition.  This  approach,  which  reduces 
concentrations  throughout  the  entire  vertical  extent  of  the  cloud  or  plume,  is  known 
as  source  depletion.  In  contrast,  OBODM  is  a  surface  depletion  model  because  it 
removes  material  just  at  the  surface  and  not  throughout  the  entire  vertical  extent  of 
the  cloud. 

OBODM  is  not  specifically  designed  to  calculate  dry  deposition,  but  can 
nevertheless  be  used  as  a  screening  tool  to  place  an  upper  bound  on  dry  deposition. 
The  worst-case  assumption  that  can  be  made  for  dry  deposition  is  that  100  percent 
of  the  pollutant  of  concern  that  is  mixed  to  the  surface  is  deposited.  The  resulting 
deposition  can  be  obtained  from  OBODM's  gravitational  deposition  option  by 
assuming  a  zero  settling  velocity  and  zero  surface  reflection  coefficient.  (For  a  gas, 
the  OBODM  user  must  assume  monodisperse  particulates  with  a  negligible  settling 
velocity  and  zero  reflection  coefficient.)  A  more  realistic  estimate  of  maximum 
possible  dry  deposition  can  be  obtained  from  Equation  (2-49).  That  is,  the  OBODM 
ground-level  dosage  predictions  can  be  multiplied  by  an  appropriate  deposition 
velocity  to  estimate  the  dry  deposition  at  each  receptor.  Because  this  approach 
effectively  assumes  that  there  are  no  dry  deposition  losses  between  the  source  and 
each  receptor  of  interest,  it  provides  the  dry  deposition  that  would  occur  if 
deposition  did  not  begin  until  each  receptor 


19 


2.6  DISPERSION  COEFFICIENTS,  TRANSPORT  WIND  SPEED,  AND  INITIAL 
SOURCE  DIMENSIONS 

2.6.1  Dispersion  Coefficients 


OBODM  uses  semi-empirical  Dugway  Proving  Ground  dispersion  coefficients 
which  directly  relate  plume  and  puff  growth  to  atmospheric  turbulence  and  vertical 
wind  shear.  These  coefficients  are  also  used  in  other  models  such  as  REEDM,  the 
FSCBG  aerial  spray  model  (Teske  et  a!.,  1993),  and  the  SHORTZ/LONGZ  complex 
terrain  dispersion  models  (Bjorklund  and  Bowers,  1982).  The  lateral  (crosswind) 
dispersion  coefficient  av  and  longitudinal  (alongwind)  dispersion  coefficient  a*  are 
each  comprised  of  atmospheric  turbulence  and  wind  shear  terms,  whereas  the 
vertical  dispersion  coefficient  az  is  a  function  of  turbulence  only. 


The  OBODM  equation  for  the  lateral  dispersion  coefficient  is  (Cramer  et  al., 
1  972) 


°y 


A  Q'x 
4.3 


1 1/2 


(2-50) 


where  oyt  is  the  turbulence  contribution  to  lateral  dispersion  and  A0’  is  the  wind- 
direction  difference  in  radians  between  the  top  and  bottom  of  the  plume  or  cloud. 
This  difference  is  normally  set  to  zero  in  hazard  analyses  to  obtain  maximum 
possible  predicted  downwind  impacts.  The  coefficient  4.3  in  Equation  (2-50)  is  the 
normalized  distance  between  the  points  on  a  Gaussian  plume  at  which  the 
concentration  is  10  percent  of  the  centerline  concentration. 


Cramer  et  al.  (1972)  derive  a  generalized  equation  for  cry,  under  the 
assumption  that  ayt  is  directly  proportional  to  the  lateral  intensity  of  turbulence  and 
varies  linearly  with  downwind  distance  x  (i.e.,  cy,  cc  x)  from  the  origin  of  an  ideal 
point  source  to  downwind  distance  xry,  beyond  which  ayt  varies  as  x“.  After 
provision  for  an  initial  source  dimension  (i.e.,  a  volume  source)  by  means  of  a 
lateral  virtual  distance,  the  resulting  ayt  equation  is 


cryt  a  ^  x,y 


x  +  xy  -  (1-a) 


a  x, 


O' 


(2-51) 
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ry 


ryJ 


(2-52) 


where 

c rA'  =  standard  deviation  of  wind  direction  angle  in  radians  [depends  on 
source  type;  see  Equation  (2-55)  below] 

xry  =  downwind  distance  from  an  ideal  point  source  over  which  rectilinear 
crosswind  .expansion  occurs  (normally  assumed  to  be  50  m) 

a  =  lateral  expansion  coefficient  (normally  assumed  to  be  0.9  for 

continuous  or  quasi-continuous  sources  and  1 .0  for  instantaneous 
sources) 

xRy  =  reference  downwind  distance  at  which  cryR  is  measured  (equals  zero  if 
o'yR  is  measured  at  the  origin) 

Note  that  the  lateral  virtual  distance  xy  given  by  Equation  (2-52)  is  not  allowed  to 
be  less  than  zero. 

The  generalized  ayt  algorithm  given  by  Equations  (2-51)  and  (2-52)  is  used  for 
both  puff  (instantaneous  source)  and  plume  (quasi-continuous  source)  dispersion, 
but  the  turbulence  inputs  and  empirical  coefficients  differ  for  the  two  source  types. 
For  quasi-continuous  sources,  the  OBODM  defaults  for  the  lateral  rectilinear 
expansion  distance  xry  and  lateral  expansion  coefficient  a  are  50  m  and  0.9, 
respectively.  If  these  default  values  are  used  with  the  SHORTZ  default  rural 
turbulence  intensities,  the  cryt  values  predicted  by  Equations  (2-51)  and  (2-52)  for  a 
point  source  closely  match  the  Pasquill-Gifford  cxy  curves  (Turner,  1970).  For 
instantaneous  sources,  the  OBODM  defaults  for  xry  and  a  are  50  m  and  1 .0, 
respectively.  However,  the  exact  value  of  xry  is  unimportant  because  Equations  (2- 
51)  and  (2-52)  reduce  to 

CTyl  =  CTa'(X  +  Xy)  (2-53) 
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OyR 
&  A' 


(2-54) 


Xy 


%  Ry 


.Thus,  lateral  puff  expansion  (neglecting  any  shear  effects)  is  normally  assumed  by 
OBODM  to  be  rectilinear  at  all  downwind  distances.  Using  the  x  =  ut  relationship, 
this  model  assumption  provides  a  good  overall  fit  to  the  available  lateral  puff 
dispersion  data  (Gifford,  1995). 


As  noted  above,  the  magnitude  of  the  cta'  used  to  calculate  CTyt  depends  on  the 
source  type.  OBODM  adjusts  cta'  from  the  measurement  time  To  to  averaging  time 
T  using  the  relationship 


oA'(T) 


o-a’(To) 


( 

vrj 


(2-55) 


where 


;  instantaneous  source  } 


;  quasi  -  continuous  source,  ts  >  tH  f 


|  MIN  (tn ,  ta  )  ;  quasi  -  continuous  source,  ts  <  tK 


(2-56) 


and  ts  is  the  cloud  stabilization  time,  tB  is  the  burn  time,  and  ta  is  the  concentration 
averaging  time.  The  stabilization  time  ts  is  obtained  by  dividing  the  distance  to  final 
rise  xr  (see  Section  2.6.3)  by  the  mean  wind  speed  u.  Note  that  T  is  not  allowed 
to  be  less  than  2.5  s.  The  t0  2  time  dependence  of  Equation  (2-55)  has  been 
suggested  by  numerous  authors  (for  example,  see  Slade,  1968;  Osipov,  1972;  and 
Wilson  and  Du,  1995)  to  account  for  the  averaging  time  variation  of  av.  See  Barr 
and  Clements  (1984)  for  a  further  discussion  of  the  rationale  for  the  use  of 
Equation  (2-55). 

The  OBODM  algorithm  for  the  vertical  dispersion  coefficient  az  is  given  by 
Equations  (2-51)  and  (2-52)  with  Ge'  (the  standard  deviation  of  the  wind  elevation 
angle  in  radians)  substituted  for  ga\  z  substituted  for  y  in  the  subscripts,  and  P 
substituted  for  a.  For  both  instantaneous  and  quasi-continuous  sources,  the  model 
defaults  for  the  vertical  rectilinear  expansion  distance  xrz  and  vertical  expansion 
coefficient  P  are  50  m  and  1 .0,  respectively.  Thus,  the  OBODM  default  values 
assume  rectilinear  expansion  analogous  to  that  of  Equations  (2-53)  and  (2-54)  at  all 
downwind  distances  for  all  stabilities.  Analyses  of  Army  field  test  data  (for 
example,  Cramer  et  al.,  1972)  suggest  that  P  can  be  approximated  as  unity  for  all 
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stabilities  for  elevated  releases,  but  depends  on  stability  for  near-surface  releases. 
Also,  it  should  be  recognized  that  the  restriction  on  mixing  at  the  top  of  the  surface 
mixing  layer  [accounted  for  by  the  multiple  reflection  terms  of  Equation  (2-28)  or 
(2-30)]  causes  the  effective  rate  of  vertical  expansion  to  decrease  with  distance  as 
the  cloud  or  plume  begins  to  fill  the  mixing  layer.  OBODM  does  not  make  any 
averaging  time  adjustments  to  a£'.  However,  when  Equation  (2-55)  is  used  to 

adjust  <ta'  to  short  averaging  times,  the  model  does  not  allow  aA'(T)  to  be  less  than 
Oe'. 


The  OBODM  algorithms  for  the  dispersion  and  deposition  of  particles  with 
appreciable  gravitational  settling  velocities  assume  that  the  particles  follow  ballistic 
trajectories  that  are  modified  by  atmospheric  turbulence.  As  the  settling  velocities 
of  particles  increase,  the  perturbations  of  their  trajectories  by  atmospheric 
turbulence  are  decreased.  OBODM  accounts  for  this  effect  by  reducing  ga'  and 
using  a  method  based  on  the  suggestions  of  Csanady  (1963,  1967).  The  effective 
turbulence  inputs  used  iq  the  avt  and  az  algorithms  for  the  nth  settling-velocity 
category  are 
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o'An  (<#)  = 


<<  (O 
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(2-57) 


<y'En  (eff)  = 
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(2-58) 


where  Pi  is  the  ratio  of  the  Lagrangian  and  Eulerian  time  scales.  Analyses  of 
Dugway  Proving  Ground  field  tests  in  which  large  drops  were  disseminated 
(Dumbauld  and  Rafferty,  1979)  have  found  that  the  Dugway  gravitational  depo¬ 
sition  algorithm  when  used  with  Equations  (2-57)  and  (2-58)  closely  reproduces  the 
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observed  deposition  patterns  if  pi  is  assumed  to  be  unity.  Although  the  OBODM 
default  for  Pi  is  unity,  other  values  may  be  entered.  If  Pi  is  entered  as  zero, 
OBODM  does  not  use  the  crossing-trajectories  correction. 

OBODM  computes  the  longitudinal  dispersion  coefficient  crx  using  the 
generalized  equation  (Dumbauld  and  Bowers,  1983) 

cr.v  =  [  cri  +  cri-  \'2  (2-59) 


<x „  =  ah  (x  +  x,) 


(2-60) 


O’.xR 

ah 


l/h 


-  XRx 


(2-61) 


(2-62) 


where 

lx 

=  longitudinal  intensity  of  turbulence,  adjusted  to  time  T  using 
Equation  (2-55) 

a,b,  and  Es 

=  empirical  constants 

CTxR 

=  standard  deviation  of  the  alongwind  concentration 
distribution  at  downwind  distance  xRx  (m) 

Au 

=  wind-speed  difference  through  the  layer  containing  the  cloud 
(m  s'1);  see  Section  2.6.2  for  the  Au  equation 

and  the  other  terms  are  as  defined  above.  (Note  that  the  longitudinal  virtual 
distance  xx  is  not  allowed  to  be  less  than  zero.) 


In  the  absence  of  user-specified  lx  values,  OBODM  uses  the  Dumbauld  and 
Bowers  (1983)  suggestion  that 

Ix  *  1.33 Jy  «  1.33 a k'  (2-63> 
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where  lv  is  the  lateral  turbulence  intensity.  The  averaging  time  adjustment  of 
Equation  (2-55)  is  applied  to  lx  with  T  equal  to  the  stabilization  time  for  a  buoyant 
source  and  2.5  s  for  a  nonbuoyant  source.  The  crossing-trajectories  correction  of 
Equation  (2-57)  is  also  applied  to  lx.  Under  the  assumption  that  the  coefficients  a 
and  b  are  both  unity,  Dumbauld  and  Bowers  (1983)  used  the  Nickoia  (1971)  data 
to  estimate  that  Es  is  approximately  0.06. 

2.6.2  Wind-Speed  Profile  and  Transport  Wind  Speed 


OBODM  assumes  that  the  height  variation  of  wind  speed  within  the  surface 
mixing  layer  can  be  described  by  the  power-law  profile 

u(z)  =  u(zR)  (2-64) 


where  ,y 

u(z)  =  wind  speed  at  height  z  above  the  surface 
u(zR)  =  wind  speed  at  height  zR  above  the  surface 

p  =  wind-profile  exponent 


The  typical  (and  default)  wind  measurement  height  zR  is  10  m. 

It  follows  from  Equation  (2-64)  that  the  wind-speed  difference  Au  required  by 
Equation  (2-62)  is 


Aw  = 


u(  Zr) 

(zn)P 


[zp2  -  zlj 


(2-65) 


where  Zi  and  Z2  are  the  heights  of  the  cloud  bottom  and  top,  respectively.  The 
bottom  of  the  cloud  zi  is  defined  as 


zi  = 


f  H O)  -  2.15  cr.  (x)  ;  H  -  2.15  ax  >2 


m 


< 


;  H  -  2.15  a.  <2  m 


(2-66) 


Similarly,  the  top  of  the  cloud  Z2  is  defined  as 


J  H(x)  +  2.15  <jz  (*)  ;  H  +  2.15  az  <  Hm 
Zl=\Hm  ;H  +  2.15  az>Hm 


(2- 
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OBODM  assumes  that  the  mean  transport  wind  speed  is  given  by 


(2-68) 


If  Equation  (2-64)  is  entered  in  Equation  (2-68),  the  resulting  expression  for  the 
transport  wind  speed  is 


-  =  u(zR)  [zVp)  -  zTr>l 

(Z2  -  Zi)  ZPR  (1  +  p) 


(2-69) 


2.6.3  Buoyant  Source  Dimensions 


Observations  of  the  buoyant  rise  of  stack  plumes  (for  example,  Briggs,  1971) 
show  that  the  visible  radius  of  the  rising  plume  can  be  approximated  by 

r(x)  =  ycAhc(x)  (2-70) 

where  Ahc  is  the  stable  or  adiabatic  plume  rise  for  a  continuous  source,  depending 
on  the  thermal  stratification.  If  this  result  is  generalized  to  the  case  of  a  source 
with  a  large  initial  horizontal  radius  rR,  Equation  (2-70)  becomes 

r(x)  =  ycAhc(x)  +  rR  (2-71) 

OBODM  sets  rR  equal  to  the  radius  of  a  circle  with  the  same  horizontal  area  as  the 
material  to  be  burned.  Under  the  assumption  of  a  bivariate  Gaussian  distribution, 
the  lateral  and  vertical  concentration  (and  dosage)  standard  deviations  at  the 
downwind  distance  xR  of  plume  stabilization  are 


CXyR 


(7-.R  = 


ycAhc(xR)  +  rR 
2.15 


(2-72) 


As  indicated  by  Equation  (2-1)  and  Equations  (2-4)  through  (2-6),  the  downwind 
distance  to  plume  stabilization  is  given  by 

_  |  nuS~m  ;  ddldz  >  0 
Xr  =  |3ix*  ;  dQldz  <  0 


(2-73) 


For  consistency  with  the  plume  rise  calculations,  xR  is  never  allowed  to  exceed 
3.5x*.  At  downwind  distances  less  than  xR,  ayR  and  crzR  are  given  by  Equation 
(2-72)  with  Ahc  evaluated  at  x  rather  than  xR. 
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Source  dimensions  are  defined  for  an  instantaneous  source  in  the  same 
manner  as  for  a  continuous  source  except  that  there  is  a  third  degree  of  freedom, 
the  alongwind  dimension.  Thus,  cloud  dimensions  at  the  downwind  distance  of 
stabilization  are  defined  by 


O'xll  ~  CTyR  ~  <7zR 


y,kh,(xR)  +  rR 

2.15 


(2-74) 


where  xR  is  given  by  Equation  (2-73)  with  xmax  [from  Equation  (2-1 1)]  substituted 
for  3.5x*  in  the  second  line.  [If  Equation  (2-7)  is  used  to  calculate  final  cloud  rise 
with  an  adiabatic  or  unstable  thermal  stratification,  xR  is  given  by  the  top  line  of 
Equation  (2-73)  with  the  stability  parameter  S  obtained  using  a  potential 

temperature  gradient  of  3.344  x  10'4  K  m'1.]  At  downwind  distances  less  than  xR, 
axR  is  set  equal  to  the  value  given  by  Equation  (2-74)  with  Ahi  evaluated  at  x  rather 

than  xR.  Using  the  approximation  that  the  average  temperature  of  the  initial 
buoyant  cloud  produced  by  a  detonation  is  1.44Ta  (Hoock  et  al.,  1987),  OBODM 
computes  rR  from 


rK 


0.89 


/  \  1/3 

3  Hc  WD 
k  4rc  Cp  pa  T„  , 


(2-75) 


where 

Hc  =  material  effective  heat  content  (cal  g'1) 

WD  =  quantity  of  material  detonated  (g) 

cp  =  specific  heat  of  air  at  constant  pressure  (0.24  cal  g'1  K'1) 
pa  =  air  density  (g  m'3) 

Ta  =  air  temperature  (K) 

If  the  release  height  h  is  entered  as  zero  for  a  detonation,  OBODM  sets  it  equal  to 
rR. 

2.6.4  Nonbuoyant  Source  Dimensions 

In  the  case  of  a  single  nonbuoyant  volume  source,  the  initial  standard 
deviations  of  the  horizontal  concentration  distributions  are  given  by 
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(2-76) 
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(2-77) 


where  W  is  the  width  of  the  short  side,  L  is  the  length  of  the  long  side,  and  b  is  the 
minimum  angle  between  the  wind  vector  and  the  long  side.  The  initial  standard 
deviation  of  the  vertical  concentration  distribution  a*o  is  given  by  the  source's  actual 
depth  or  vertical  dimension  divided  by  2.15  for  a  surface-based  (h  =  0)  source  and 
by  4.3  for  an  elevated  (h>0)  source. 

As  noted  in  Section  1,  OBODM  represents  a  line  source  by  a  series  of 
overlapping  volume  sources.  The  number  of  volume  sources  used  to  represent  a 
line  source  is  given  by  the  line's  length  L  divided  by  its  width  W,  rounded  up  to  the 
nearest  integer.  For  each  volume  source  used  to  represent  a  line  source,  axo  is 
given  by  Equation  (2-76)  with  L  replaced  by  W,  while  ay°  is  given  by  Equation 
(2-77)  with  L  replaced  by  W  and  4.3  replaced  by  2.3548.  The  initial  standard 
deviation  of  the  vertical  concentration  distribution  is  defined  in  the  same  manner  as 
for  a  volume  source. 

2.7  COMPLEX  TERRAIN  SCREENING  PROCEDURES 

OBODM  uses  complex  terrain  screening  procedures  that  are  based  on  the 
procedures  used  by  the  SHORTZ/LONGZ  complex  terrain  dispersion  models 
(Bjorklund  and  Bowers,  1982).  Although  these  procedures  have  been  used  with 
some  success  in  modeling  industrial  stack  emissions  in  complex  terrain  (for 
example,  Cramer  et  al.,  1975),  they  are  likely  to  overestimate  the  ground-level 
impacts  of  OB/OD  releases  in  many  situations.  The  OBODM  complex  terrain 
procedures  also  do  not  address  the  buildup  of  OB/OD  emissions  that  may  become 
trapped  within  deep  valleys  during  air  stagnation  episodes. 

When  OBODM  is  applied  in  complex  terrain,  the  receptor  height  z  above 
ground  level  must  be  zero,  and  the  model  cannot  be  used  to  calculate  concen¬ 
tration,  dosage,  or  deposition  for  particulates  with  appreciable  settling  velocities 
because  the  assumptions  upon  which  Equations  (2-30),  (2-43),  (2-44),  and  (2-46) 
through  (2-48)  are  based  are  violated.  The  OBODM  complex  terrain  procedures 
modify  the  flat  terrain  model  methodologies  described  in  the  preceding  sections  by 
defining  effective  cloud  or  plume  heights  and  mixing  depths.  The  following 
assumptions  are  made  in  complex  terrain: 
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a.  The  actual  top  of  the  surface  mixing  layer  extends  over  the  calculation 
grid  at  a  constant  height  above  mean  sea  level.  The  actual  top  of  the 
surface  mixing  layer  should  not  be  confused  with  the  effective  top  of  the 
surface  mixing  layer,  which  is  a  mathematical  device  used  to  preclude 
violations  of  the  Second  Law  of  Thermodynamics  when  clouds  or  plumes 
pass  over  elevated  terrain. 

b.  The  centroid  of  a  cloud  or  axis  of  a  plume  contained  within  the 
surface  mixing  layer  remains  at  the  stabilization  height  above  mean 
sea  level,  and  the  cloud  or  plume  may  impact  elevated  terrain  under 
stable,  neutral,  or  unstable  conditions. 

c.  Clouds  or  plumes  that  stabilize  above  the  top  of  the  surface  mixing  layer 

do  not  contribute  to  significant  ground-level  concentrations  at  any 
receptor,  including  receptors  that  are  above  the  top  of  the  surface  mixing 
layer.  ,-/ 

The  OBODM  transitional  (distance-dependent)  plume/cloud  rise  algorithms  are 
not  used  in  complex  terrain.  That  is,  the  final  plume/cloud  rise  is  assumed  to  apply 
at  all  downwind  distances,  including  distances  less  than  the  distance  to 
plume/cloud  stabilization.  The  height  Ho  above  mean  sea  level  of  the  axis  of  a 
stabilized  plume  or  centroid  of  a  stabilized  cloud  is  given  by  the  sum  of  the  height  H 
above  ground  level  of  the  stabilized  plume  or  cloud  and  the  elevation  zs  above  mean 
sea  level  of  the  OB/OD  source. 

In  order  to  determine  the  fraction  of  the  plume  or  cloud  that  is  contained 
within  the  surface  mixing  layer,  the  mixing  depth  Hms  at  the  source  is  substituted 
for  Hm  in  Equation  (2-14).  The  mixing  depth  at  the  source  is  given  by 

H  —  H  +  z  -  z  (2-78) 

7  i  ni.s  lima  '  &  a  Av 


where  Hma  is  the  mixing  depth  measured  at  a  point  with  elevation  za  above  mean 
sea  level.  Note  that  Hms  is  also  used  in  Equation  (2-16)  to  determine  the  effective 
plume  or  cloud  rise  for  the  material  in  the  surface  mixing  layer. 

The  effective  height  of  the  cloud  centroid  or  plume  axis  above  a  receptor  with 
elevation  zr  above  mean  sea  level  is  defined  as 


H'(zr) 


Ho  -  zr  ;  Ho  -  zr  >ol 
0  ;  Ho  -Zr  <0 J 


(2-79) 


The  effective  height  H'(zr)  is  substituted  for  H  in  the  vertical  term  of  Equation  (2- 
30).  The  effective  mixing  depth  above  the  same  receptor  is  given  by 
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The  SHORTZ/LONGZ  models  allow  the  user  to  assume  that  the  wind  speed  is 
a  function  of  height  above  the  ground  surface  or  height  above  mean  sea  level. 
Rather  than  allow  this  user  choice,  OBODM  always  assumes  that  the  wind  speed  is 
a  function  of  the  height  above  the  local  ground  level.  The  wind  profile  and 
transport  wind  speed  equations  in  Section  2.6.2  are  used  without  modification 
except  that  Hma  is  substituted  for  Hm  in  the  definition  of  Z2  in  Equation  (2-67). 


SECTION  3.  DEFAULT  METEOROLOGICAL  INPUTS 


The  minimum  OBODM  meteorological  inputs  consist  of  the  wind  speed  and 
wind  direction  at  10  m  above  the  surface,  the  air  temperature,  and  either  the 
Pasquill  stability  category  or  the  Net  Radiation  Index  (NRI).  With  this  minimum 
information,  the  program  will  provide  default  values  for  all  other  meteorological 
inputs.  However,  the  user  is  encouraged  to  use  measured  or  site-specific 
climatological  values  whenever  available. 

Table  2  summarizes  the  meteorological  conditions  represented  by  the  Pasquill 
stability  categories.  As  shown  by  Table  3,  Turner  (1964)  defines  each  Pasquill 
stability  category  by  the  combination  of  the  10-m  wind  speed  and  NRI.  Table  4 
lists  the  steps  used  to  estimate  the  NRI  from  the  solar  altitude  angle  and  standard 
surface  observations  of  cloud  cover  and  ceiling  height.  The  OBODM  default 
meteorological  inputs  are  assigned  based  on  the  combination  of  the  10-m  wind 
speed  and  NRI.  If  the. Pasquill  stability  category  rather  than  the  NRI  (or  the 
information  needed  to  compute  the  NRI)  is  entered,  OBODM  uses  Table  5  to 
estimate  the  NRI. 


Table  2.  Meteorological  Conditions  Represented  by  the  Pasquill  Stability 
Categories. 


Pasquill  Stability  Category 

Meteorological  Conditions 

A 

Very  Unstable 

B 

Unstable 

C 

Slightly  Unstable 

D 

Neutral 

E 

Stable 

F 

Very  Stable 
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Table  3.  Pasquill  Stability  Categories  as  a  Function  of  Wind  Speed  and 
Net  Radiation  Index  (Turner,  1964). 


Wind  Speed 

Net  Radiation  Index 

at  1 0  m  3 
(kt) 

4 

3 

2 

1 

0 

-1 

-2 

0,1 

A 

A 

B 

C 

D 

E 

F 

2,3 

A 

B 

B 

C 

D 

F 

F 

4,5 

A 

B 

C 

D 

D 

E 

F 

6 

B 

B 

C 

D 

D 

E 

F 

7 

B  '  ' 

U 

B 

B 

C 

D 

D 

D 

E 

8,9 

C 

C 

D 

D 

D 

E 

10 

C 

C 

D 

D 

D 

D 

E 

1 1 

\ 

C 

C 

D 

D 

D 

D 

D 

112 

c 

D 

D 

D 

D 

D 

D 

3  The  Turner  (1964)  scheme  uses  wind  speeds  in  knots  (1  kt  =  0.514  m/s). 
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Table  4.  Procedures  Used  to  Calculate  Net  Radiation  Indices  from  Cloud  Cover 
and  Ceiling  Observations  and  the  Solar  Altitude. 


1 .  If  the  total  cloud  cover  is  10/10  (overcast)  and  the  cloud  ceiling  is  less  than 
21  34  m  (7000  ft),  the  net  radiation  index  (NRI)  is  0  (whether  day  or  night). 

2.  For  night  (between  sunset  and  sunrise): 

(a)  If  the  total  cloud  cover  is  less  than  or  equal  to  4/10,  use  NRI  equal  to  -2. 

(b)  If  the  total  cloud  cover  is  greater  than  4/10,  use  NRI  equal  to  -1 . 

3.  For  daytime: 

(a)  Determine  the  insolation  class  number  from  the  solar  altitude 
(a)  as  follows:  ,j 

SOLAR  ALTITUDE  INSOLATION  CLASS 
(deg)  ,  NUMBER 

a  >  60  4 

60  -  a  >  35  3 

35  -  a  >  15  2 

1  5  -  a  1 

(b)  If  the  total  cloud  cover  is  less  than  or  equal  to  5/10,  the  NRI  is 
equal  to  the  insolation  class  number. 

(c)  If  the  total  cloud  cover  is  greater  than  5/10,  the  NRI  is  equal  to  the 
insolation  class  number  modified  using  the  following  five  steps: 

(i)  If  the  ceiling  is  less  than  2134  m  (7000  ft),  subtract  2. 

(ii)  If  the  ceiling  is  greater  than  or  equal  to  21 34  m  (7000  ft)  but  less  than 
4877  m  (16,000  ft),  subtract  1. 

(iii)  If  the  total  cloud  cover  equals  10/10,  subtract  1.  (This  will  only  apply 
to  ceilings  greater  than  2134  m  (7000  ft)  since  cases  with  10/10  coverage 
below  2134  m  are  considered  in  item  1  above.) 

(iv)  If  the  insolation  class  number  has  not  been  modified  by  steps  (I),  (ii), 
or  (iii)  above,  assume  the  NRI  is  equal  to  the  insolation  class  number. 

(v)  If  the  modified  insolation  class  number  is  less  than  1,  let  the  NRI 
equal  1 . 
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Table  5.  Net  Radiation  Index  as  a  Function  of  Pasquill  Stability  Category  and  Wind 


Speed. 


Pasquill 

Stability 

Category 

10-m  wind  Speed  (kt) 

0-1 

>1-3 

>3-5 

'  >5-6 

>6-7 

>7-9 

>9-10 

>10-11 

>11 

A 

4 

4 

4 

4 

4 

4 

4 

4 

4 

B 

2 

3 

3 

3 

3 

4 

4 

4 

4 

C 

1 

1 

2 

2 

2 

3 

4 

4 

4 

D 

0 

0 

1 

1 

0 

0 

0 

0 

0 

E 

-2 

-2 

-1 

-1 

0 

0 

0 

0 

0 

F 

-2 

-2 

-2 

-2 

-1 

-1 

-1 

0 

0 

it 


Table  6.  Default  Standard  Deviations  of  the  Wind  Direction  Angle  in  Degrees 
Measured  Over  600  s  at  a  Height  of  10  m  Above  Open  Terrain 
(Dumbauld  and  Bowers  1983). 

Net  Radiation  Index 

Wind  Speed 
(m/s) 

4 

3 

2 

1 

0 

1,-2 

0  -  1 

26.0 

22.0 

21.0 

15.0 

1 1.0 

4.0 

>1-3 

22.0 

20.0 

16.0 

1 1.0 

10.0 

4.0 

>3-5 

19.0 

16.0 

13.0 

1 1.0 

9.0 

7.0 

>5-7 

15.0 

13.0 

12.0 

9.0 

8.0 

6.7 

>7 

13.0 

12.0 

10.0 

8.0 

7.0 

6.7 

Table  7.  Default  Standard  Deviations  of  the  Wind  Elevation  Angle  in  Degrees 
at  10  m  Above  Open  Terrain  (Dumbauld  and  Bowers,  1983). 

1 0-m  Wind  Speed 

(m/s)  4 

Net  Radiation  Index 

3  2  10 

-1,-2 

0  -  1 

8.7 

8.7 

7.0 

5.0 

4.0 

3.0 

>1-3 

7.0 

7.0 

6.0 

4.0 

3.5 

3.0 

>3  -  5 

7.0 

6.0 

5.5 

5.3 

5.0 

5.0 

>5  -  7 

6.0 

5.3 

5.0 

4.8 

4.8 

4.7 

>7 

5.3 

5.0 

4.5 

4.5 

4.5 

4.5 
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The  OBODM  turbulence  inputs  consist  of  the  standard  deviations  of  the  wind 
direction  and  elevation  angles  and  the  longitudinal  turbulence  intensity.  (The 
standard  deviations  of  the  wind  direction  and  elevation  angles  in  radians  are 
respectively  equivalent  to  the  lateral  and  vertical  turbulence  intensities.)  The 
default  values  for  the  standard  deviations  of  the  wind  direction  and  elevation  angles 
are  listed  in  Tables  6  and  7,  respectively.  These  values,  which  are  intended  to  be 
representative  of  measurements  made  at  10  m  above  a  surface  with  a  roughness 
length  zo  of  about  10  cm,  were  adapted  by  Dumbauld  and  Bowers  (1983)  from  the 
values  suggested  by  Dumbauld  (1982).  The  Dumbauld  (1982)  turbulence 
intensities,  which  are  for  a  measurement  height  of  5  m  above  open  terrain,  were 
based  on  analyses  of  turbulence  measurements  at  White  Sands  Missile  Range,  New 
Mexico  (Swanson  and  Cramer,  1965);  Round  Hill,  Massachusetts  (Cramer  et  al., 
1966);  and  Dugway  Proving  Ground,  Utah  (Cramer  et  al.,  1972).  Dumbauld  (1982) 
also  considered  the  turbulence  profiles  shown  by  Lumley  and  Panofsky  (1964)  in 
arriving  at  his  recommehded  values.  The  OBODM  default  longitudinal  turbulence 
intensities  are  1 .33  times  the  wind  direction  standard  deviations  in  Table  6 
(Counihan,  1  975).  If  the  user  enters  the  standard  deviation  of  the  wind  direction, 
but  not  the  longitudinal  turbulence  intensity,  OBODM  sets  the  longitudinal 
turbulence  intensity  equal  to  1 .33  times  the  user-specified  wind  direction  standard 
deviation. 

OBODM  does  not  make  any  height  adjustments  to  the  default  turbulence 
inputs  in  Tables  6  and  7.  However,  the  program  makes  surface  roughness 
adjustments  if  the  user  enters  a  nonzero  roughness  length  zo.  The  roughness 
adjustment  is 

aA(z0)  =  a A  (Table  6) 


cr i:  ( z„)  =  (7 1;  (Table  7) 


where  zo  is  in  centimeters  and  aA  and  °e  are  the  standard  deviations  of  the  wind 
direction  (azimuth)  and  elevation  angles,  respectively.  Equations  (3-1)  and  (3-2)  are 
taken  from  Hanna  et  al.  (1977). 

Table  8  gives  the  OBODM  default  wind-profile  (power  law)  exponents.  These 
exponents  were  adapted  by  Dumbauld  and  Bowers  (1983)  from  the  coefficients 
suggested  by  Dumbauld  (1982),  which  were  based  on  analyses  of  wind  speed 
profiles  at  White  Sands  Missile  Range  and  Dugway  Proving  Ground. 


10 


0.2 


Zu_ 

10  J 


\  02 


(3-1) 


(3-2) 
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Table  8.  Default  Wind  Speed  Power  Law  Exponents  Over  Open  Terrain 
(Dumbauld  and  Bowers,  1983). 


Net  Radiation  Index 

10-m  Wind  Speed 
(m/s) 

4 

3 

2  1  0 

-1,-2 

0  -  1 

0.10 

0.10 

0.20 

0.20 

0.25 

0.40 

>1-3 

0.10 

0.10 

0.20 

0.20 

0.25 

0.30 

>3  -  5 

0.10 

0.10 

.  0.15 

0.15 

0.20 

0.25 

>5  -  7 

0.10 

0.10 

0.10 

0.10 

0.10 

0.20 

>7 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

„•/ 


Table  9.  Default  Vertical  Potential  Temperature  Gradients  in  Degrees  Kelvin 
Per  Meter  (Bjorklund  and  Bowers,  1982). 


_ Pasquill  Stability  Category _ 

10-m  Wind  Speed 

(m/s)  A  B  C  D  E  F 


Relative  Humidity  >70% 


0  -  1 

0.000 

0.000 

0.000 

0.015 

0.030 

0.035 

>1-3 

0.000 

0.000 

0.000 

0.010 

0.020 

0.025 

>3  -  5 

0.000 

0.000 

0.000 

0.005 

0.015 

0.015 

>5  -  7 

0.000 

0.000 

0.000 

.0.003 

0.010 

0.010 

>7 

0.000 

0.000 

0.000 

0.003 

0.003 

0.003 

Relative  Humidity  - 

70% 

0  -  1 

0.000 

0.000 

0.000 

0.020 

0.030 

0.040 

>1-3 

0.000 

0.000 

0.000 

0.010 

0.020 

0.030 

>3  -  5 

0.000 

0.000 

0.000 

0.005 

0.010 

0.020 

>5  -  7 

0.000 

0.000 

0.000 

0.000 

0.005 

0.010 

>7 

0.000 

0.000 

0.000 

0.000 

0.000 

0.005 
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Table  10.  Default  Surface  Mixing  Layer  Depths  in  Meters  for  Open  Terrain 
(Dumbauld  and  Bowers,  1983). 


1 0-m  Wind  Speed 
(m/s) 

Net  Radiation  Index 

4 

3 

2  1 

0 

-1,-2 

0  -  1 

2500 

2000 

1000 

500 

100 

30 

>1-3 

2200 

1800 

1200 

600 

200 

100 

>3-5 

1800 

1500 

1200 

600 

300 

200 

>5  -  7 

1500 

1200 

1000 

600 

300 

200 

>7 

1200 

1000 

700 

500 

300 

200 

ij 

The  OBODM  default  vertical  potential  temperature  gradients  and  mixing  layer 
depth  are  listed  in  Tables  9  and  10,  respectively.  The  potential  temperature 
gradients  in  Table  9  are  the  default  inputs  used  by  the  SHORTZ  and  LONGZ 
diffusion  models  (Bjorklund  and  Bowers,  1982).  The  potential  temperature 
gradients  for  high  relative  humidities  are  principally  taken  from  a  modeling  study  of 
sulfur  dioxide  (SOz)  air  quality  in  Allegheny  County,  Pennsylvania  (Cramer  et  al., 
1975).  The  potential  temperature  gradients  for  the  lower  relative  humidities  are 
principally  based  on  measurements  made  at  Dugway  Proving  Ground.  The  mixing 
depths  in  Table  10,  which  are  taken  from  Dumbauld  and  Bowers  (1983),  are 
generally  representative  of  noncoastal  regions  in  the  western  United  States  (for 
example,  see  Holzworth,  1972). 
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SECTION  4.  EMISSIONS  FACTOR  DATABASE 


The  OBODFUEL.OBD  file  contains  a  database  of  source  information  required 
by  OBODM  for  36  propellant,  explosive,  and  pyrotechnic  energetic  fuels/  materials. 
For  each  fuel/material,  the  file  contains  a  list  of  gases,  volatile  organic  compounds, 
semivolatile  organic  compounds,  and  metals  produced  by  open  burning  or  open 
detonation.  For  each  gas,  compound,  or  metal,  the  file  provides  the  molecular 
weight,  density  at  1  atmosphere  and  20  °C,  and  emissions  factor.  As  discussed 
below,  the  emissions  factors  were  derived  from  OBOD  experiments  conducted  at 
U.S.  Army  Dugway  Proving  Ground.  The  OBODM  user  can  modify  the 
OBODFUEL.OBD  database  file  by  entering  updated  emissions  data  or  emissions 
data  for  new  fuels  or  materials. 

The  emissions  factors  in  the  OBODFUEL.OBD  file  were  developed  from  the 
emissions  data  in  the  1997  OBOD  Database  Lite  VI  .2  program.  The  majority  of 
the  information  in  the  1997  Database  Lite  VI. 2  comes  from  experiments  in 
Dugway's  BangBox™  facility  in  which  small  amounts  of  the  fuel  or  material  were 
burned  or  detonated,  and  the  resulting  gases,  compounds,  and  metals  were 
sampled  and  assayed  to  quantify  the  emissions.  Limited  data  were  also  obtained 
from  aircraft  sampling  of  open  air  OBOD  events  during  earlier  OBOD  programs  at 
Dugway.  Depending  on  the  fuel/material  and  combustion  product,  the  1997  OBOD 
Database  Lite  VI  .2  contains  emissions  factors  for  one  to  three  trials.  If  emissions 
factors  are  available  for  only  one  trial,  those  emissions  factors  are  provided  in  the 
OBODFUEL.OBD  file.  If  emissions  factors  are  available  for  two  trials,  OBODFUEL. 
OBD  contains  the  highest  emissions  factor  for  each  gas,  compound,  or  metallic 
element  as  a  safe-sided,  "worst-case"  estimate.  If  emissions  factors  are  available 
for  three  trials,  OBODFUEL.OBD  contains  the  average  emission  factor  for  each 
combustion  product.  In  some  cases,  two  different  sampling  techniques  were  used 
for  the  same  constituent.  In  these  cases,  OBODFUEL.OBD  provides  the  highest  of 
the  average  emissions  factors  obtained  for  that  constituent  using  the  two  sampling 
methods.  If  a  gas,  compound,  or  metal  was  below  the  detection  limit  for  all  trials, 
that  constituent  is  not  included  in  the  model  database  file. 

The  format  of  the  OBODFUEL.OBD  file  is  the  MS-DOS  text  mode.  The  first 
line  for  each  of  the  36  propellant,  explosive,  and  pyrotechnic  energetic  fuels  and 
materials  contains  the  name  of  the  fuel/material,  the  estimated  heat  content  and 
burn  rate,  and  the  number  (N)  of  combustion  products  for  which  emissions  data  are 
provided.  This  line  is  followed  by  N  lines,  each  listing  the  molecular  weight, 
density,  and  emissions  factors  for  one  of  the  N  products  of  combustion.  Thus, 
OBODFUEL.OBD  contains  N  +  1  lines  per  propellant,  explosive,  or  pyrotechnic 
energetic  fuel  or  material.  The  effective  heat  contents  and  burn  rates  are  not 
known  for  many  materials.  For  those  fuel/materials  with  unknown  heat  contents, 
OBODFUEL.OBD  assigns  a  value  of  1000  cal/g.  This  value  should  underestimate 
the  actual  heat  release  and  buoyant  rise  for  most  fuels/materials,  resulting  in  a 
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tendency  toward  overestimation  of  maximum  ground-level  impacts.  No  burn  rates 
are  provided  if  the  burn  rate  is  not  known.  Also,  in  some  cases  the  molecular 
weight  and/or  density  could  not  be  found  and  are  not  provided. 
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SECTION  5.  APPENDICES 


APPENDIX  A.  DERIVATION  OF  EQUATION  (2-1 1 ) 


Briggs  (1970)  assumes  that  only  eddies  with  sizes  on  the  same  order  as  the 
plume  radius  are  effective  in  mixing  ambient  air  into  the  plume.  These  eddies  are  in 
the  inertial  subrange,  which  can  be  characterized  by  the  eddy  dissipation  rate  e. 

For  eddies  on  the  order  of  the  plume  radius  r,  Briggs  assumes  that 

—  =Bsl/3r,/3  (A-1) 

dt  r 


where  p  is  a  dimensionless  constant. 

i/ 

The  closure  assumption  used  by  Briggs  is 

r  =  yz 


(A-2) 


where  y  is  the  entrainment  coefficient  and  z  is  the  plume  height,  which  is  given  by 
Briggs'  equation  for  the  rise  of  a  buoyant  plume  in  an  adiabatic  atmospheric.  It 
follows  from  Equation  (A-2)  that 


dr 

dt 


dz 

=  yjt=r» 


(A-3) 


where  w  is  the  plume's  vertical  velocity.  Briggs  combines  Equations  (A-1)  and  (A- 
3)  and  solves  for  the  distance  x*  at  which 


dr 

dt 


yw  =  Pe1'3  rl/3 


(A-4) 


This  distance  is  given  by 


x* 


3/3 


.1/3  \-9/3 


(PsU3) 


(A-5) 


where  F  is  the  plume’s  buoyancy  flux.  Briggs  then  reasons  that  maximum  or  final 
plume  rise  is  attained  at  3.5  x*. 


A-1 


Equation  (A-4)  also  forms  the  basis  for  estimating  the  distance  to  final  rise  for  a 
buoyant  cloud  in  an  adiabatic  or  unstable  atmosphere.  Neglecting  the  initial 
radius  rR,  Dumbauld  et  al.  (1973)  give  the  cloud  height  at  downwind  distance  x  in 
an  adiabatic  atmosphere  as 


Z,= 


2  F,  x2 

3  T 

r,  « 


-1 1/4 


(A-6) 


where  Fi  is  the  cloud's  buoyancy  flux  and  u  is  the  mean  wind  speed.  The  closure 
assumption  of  Equation  (A-2)  then  gives 


r  =  r,z,  = 


2F,  y,  x2 


“1 1/4 


The  vertical  velocity  w  and  wind  speed  u  are  related  by 


dx  dz ,  dz, 

w=  — - —  =  u  — 

dt  dx  dx 


Combination  of  Equations  (A-6)  and  (A-8)  yields 

F,  u2 


w  = 


-|  1/4 


(2  r,y 


X 


-1/2 


Consequently, 


r,  w  = 


f,  r,  u2 

8x2 


“1 1/4 


(A-7) 


(A-8) 


(A-9) 


(A-10) 


If  Equations  (A-7)  and  (A-10)  are  substituted  into  Equation  (A-4),  the  resulting 
distance  is 


X,  =  [/?*"]■"  (T)  (F,r,r«  IA-11) 

Use  of  Equation  (A-1 1)  to  determine  the  downwind  distance  of  final  cloud  rise 
3.5x*  requires  estimates  of  P  and  Y.  Following  Briggs  (1970),  P  is  conservatively 
estimated  as  unity  and  8  is  assumed  to  be  given  by 


A-2 


£  = 


0.068  u  /  z  ;  z  <  100  m 
0.00068  u  ;  z  >  100  m 


(A-12) 


where  z  is  height  above  the  ground.  For  yi  equal  to  0.64  and  z  equal  to  Zi  and  less 
than  or  equal  to  100  m,  substitution  of  Equation  (A-12)  into  Equation  (A-1 1)  gives 
the  distance  to  final  cloud  rise  as 


3.5x*  =  5z)/2  F; 


r*  1/4  .1/2 


u 


(A-1 3) 


If  Equations  (A-6)  and  (A-1 3)  are  combined  and  solved  for  Zi,  the  predicted  final  rise 
is 


z,  (3.5x*)  = 


3/2 


50  F, 

V ]  » 


1/3 


(A-1 4) 


Comparison  of  Equations  (A-6)  and  (A-1 4)  shows  that  Equations  (A-1 4)  can  be 
rewritten  as 


where 


*/  (Xmax)  = 


IF, 


r, 


1/4 


(A-1 5) 


xmax  =12  F 


1/2  um 


(A-1 6) 


As  indicated  by  Equation  (A-12),  Equations  (A-1  5)  and  (A-1 6)  are  not  valid  for 
cloud  rises  above  100  m.  If  the  left-hand  side  of  Equation  (A-1 4)  is  set  equal  to 
100  m  and  yi  is  set  equal  to  0.64,  the  resulting  range  of  applicability  of  Equations 
(A-1  5)  and  (A-1  6)  is 


F,  <  300  u2n 


(A-1 7) 


For  larger  values  of  Fi,  the  second  line  of  Equation  (A-12)  should  be  used  in 
Equation  (A-1 1).  The  resulting  distance  to  final  cloud  rise  is 

3.5  x*  =50  F‘l4  uU2  (A-1 8) 


The  generalized  cloud  rise  equations  under  adiabatic  or  unstable  conditions  then 
become 


A- 3 


Ah  = 


(A- 19) 


2  Fi  x2 
Y3,  u- 


-1//V 


;  X^Xmnx  | 

;  X  >  Xmax  J 


(A-20) 


12  F]n  um  ;  F,  <  300  um  1 
50  F//4  uV2  ;  F,  >  300  w2/3  J 


(A-21) 
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